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Abstract. Is there a connection between the branch point singularity at the particle 
emission threshold and the appearance of cluster states which reveal the structure of 
a corresponding reaction channel? Which nuclear states are most impacted by the 
coupling to the scattering continuum? What should be the most important steps in 
developing the theory that will truly unify nuclear structure and nuclear reactions? 
The common denominator of these questions is the continuum shell-model approach 
to bound and unbound nuclear states, nuclear decays, and reactions. 
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1. Introduction 

The nuclear shell model was proposed almost sixty years ago [U [2] . Soon afterwards, the 
interacting shell model (ISM) was developed by Lane [3], Kurath [I], and others (see [5] 
for a recent review). The ISM-based description of an evolution of the nucleonic coupling 
scheme from the LS to jj coupling with increasing mass number [31 El E], provided 
foundations of modern nuclear structure theory and helped to understand and categorize 
a wealth of data on nuclear levels, moments, collective excitations, electromagnetic and 
j3 decays, and various particle decays [7]. 

In its traditional form, ISM describes the nucleus as a closed quantum system: 
nucleons occupy bound, hence well localized, single-particle orbits of an infinite (e.g., 
harmonic oscillator) potential and are isolated from the environment of unbound 
scattering states that are not square integrable. Since the scattering continuum is 
not considered explicitly, the presence of branch points (decay thresholds) and double- 
poles of the scattering matrix (S'-matrix) is neglected. The divide between the discrete 
states and the scattering continuum, i.e., the focus on one or another, has unfortunately 
become a kind of paradigm. In the long term, this has led to an artificial separation 
of nuclear structure from nuclear reactions, and hindered a deeper understanding of 
nuclear properties. Indeed, many structural properties of the nucleus are determined 
by means of nuclear collisions. Hence, the knowledge of nuclear structure depends on 
nuclear reactions and vice versa, and this cries out for a unified theoretical framework. 

The first attempt in this direction came from Feshbach [HI [9] who expressed the 
collision matrix of the optical model using matrix elements of the Hamiltonian. This 
development gave a strong push to the ISM approach to nuclear reactions [HE [TTJ [121 [13] 
(see introduction in Ref. [H] for a detailed historical account) and led to various 
formulations of the continuum shell model (CSM) [151 EE El HI HI 120] • A modern 
version of CSM in the Hilbert space, the Shell Model Embedded in the Continuum 
(SMEC) [HI [21] , provides a unified description of the structure and reactions with 
up to two nucleons in the scattering continuum using realistic ISM Hamiltonians. 
Nevertheless, the fully symmetric treatment of bound and scattering states in the 
multiparticle wave function is still too ambitious a goal |22j. 

A different attempt to formulate the ISM for open quantum systems has been 
proposed recently [231 1^1 EE 126] within the Berggren ensemble [27]. The resulting 
complex-energy open quantum system extension of the ISM, the Gamow Shell Model 
(GSM), can be conveniently formulated in the Rigged Hilbert Space (Gel'fand triple) 
[281 [29] . which encompasses Gamow states [30l [31], and is suitable for extending the 
quantum mechanics into the domain of time-asymmetric processes (e.g. decays). The 
GSM offers a fully symmetric treatment of bound, resonance, and scattering single- 
particle states but, until now, has been primarily used in the context of nuclear structure. 
(For a recent review of the complex-energy shell model, see Ref. |32j.) 

In this paper, we shall draw attention to salient threshold effects related to generic 
properties of branch points and ^-matrix poles that are expected to impact structural 
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properties of nuclei. In Sec. 12. 2\ we shall use the phase rigidity indicator to demonstrate 
effects of the branch point at the particle emission threshold and the exceptional 
point (the double-pole of the 5-matrix) on the configuration mixing in the continuum. 
Section 12.31 discusses the stability of ISM eigenfunctions in the neighborhood of the 
reaction threshold. Finally, in Sec. 12.41 we outline necessary developments to achieve 
a unified theory of structure and reactions offering a symmetric treatment of bound, 
resonance, and scattering states. 

2. Salient near-threshold phenomena 

At low excitation energies, well-bound nuclei can be considered as closed quantum 
systems, well described by the standard ISM or its modern versions such as the no-core 
shell model (33j |34j [35]. Moving towards drip lines, or higher in excitation energy, the 
continuum coupling becomes gradually more important, changing the nature of weakly 
bound states. (Properties of unbound states are impacted by couplings to reaction 
channels.) In this regime, the chemical potential has a magnitude similar to the pairing 
gap; hence, the system is dominated by many-nucleon correlations which no longer 
can be considered as small perturbations atop the average potential [36]. Many-body 
states in neighboring nuclear systems with different proton and neutron numbers become 
interconnected via continuum, forming correlated domains (clusters) of quantum states. 
At present, little is known about the basic features of such clusters, e.g. a typical 'cluster 
size' and its dependence on binding properties of participating states. 

In standard network problems, like the percolation problem [371 EE] , different phases 
of correlated domains are formed depending on the mean network activity with respect to 
some threshold. Analogously, one expects that the key to understanding the formation of 
correlated domains of quantum states is a near-threshold behavior of a single network 
member. In the correlation-dominated regime of a nuclear system, the percolation 
threshold(s) should be identified with reaction threshold(s). 

What can be said about properties of many-body systems in the narrow range 
of energies around the reaction threshold? Are those properties independent of any 
particular realization of the Hamiltonian? Below, we shall try to address these questions 
in the framework of SMEC, but essentials of this discussion are valid in the GSM 
approach as well. 

2.1. The main ingredients of the continuum shell model 

To facilitate further discussion, we shall now recall fundamentals of SMEC. The 
Fock space of an A— particle system consists of the set of square-integrable functions 
Q = {$S A) }, internal space described by the standard ISM, and the space of embedding 
scattering states V = {CI;}- These two orthogonal and complementary subspaces of 
the Fock space are obtained by the projection formalism. In one-particle continuum 
problems, reaction channels 'c' correspond to the motion of an unbound nucleon in a 
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Figure 1. At low excitation energy, well-bound nuclei close to the stability valley can 
be described as closed quantum systems. When moving towards the particle drip lines, 
nuclei with different proton and neutron numbers become interconnected via continuum 
space, i.e., reaction channels, forming correlated domains of quantum states. 



state (Ij) relative to the daughter nucleus in some ISM eigenstate $Jf An open 
quantum system description of 'internal' dynamics, i.e., in Q, includes couplings to 
the 'environment' of scattering states and decay channels and is given by the energy- 
dependent effective Hamiltonian: 

H QQ (E) = H QQ + W(E), (1) 

where Hqq is the ISM Hamiltonian (the closed quantum system Hamiltonian) in the 
internal space Q and 

W(E) = H QP G { +\E)H PQ (2) 

is the energy- dependent term describing both virtual particle excitations and irreversible 
decays to the environment of reaction channels. In the latter expression, Gp~\E) 
is the Green's function for the motion of a single nucleon in the continuum, and 
E is the nucleon's energy (i.e., the scattering energy). The external mixing of two 
ISM states i and j due to W(E) consists of the Hermitian principal value integral 
W R (E) describing virtual excitations to the continuum and the anti-Hermitian residuum 
W I (E) = — (z/2)VV T in a dyadic form that represents decay out of the internal space 
Q. The M x A matrix V = {Vf} denotes the amplitudes connecting the ISM state $j 
(i — 1, • • • , M) to the reaction channel c (c = 1, ■ • • , A) [22]. 

At energies below the lowest reaction threshold, the effective Hamiltonian is 
Hermitian, i.e., W 1 = 0. Above the first threshold, the non-Hermitian part of Hqq(E) 
describes the irreversible decay from the internal space. The effective Hamiltonian in 
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this case becomes complex symmetric. Each eigenstate of Hqq(E) is coupled to states 
in neighboring nuclei via a network of reaction channels, either closed or open. The 
contribution of different reaction channels to the total continuum coupling is highly 
non-uniform and spans over a considerable range of excitation energies [39]. Since the 
Hamiltonian ([Q) depends explicitly on the energy, it is highly non-linear. Moreover, 
the continuum-coupling term generates effective many-body interactions in the internal 
space, even if it is originally two-body in the full space. 

Effects of the resulting many-body interactions have been studied by considering 
the continuum-coupling correction to the binding energy for the chains of oxygen and 
fluorine isotopes j39j HO]. It has been found that the induced many-body interactions 
explain a significant shift in the neutron drip line for fluorine isotopes as compared to the 
oxygens. Another explanation, explicitly invoking effective three-body interaction, has 
recently been suggested in Ref. jH]. These two scenarios could be difficult to distinguish. 

Eigenstates of the open quantum system Hamiltonian 'Hqq(E) are biorthogonal. 
The left \^fj) and right eigenstates have complex conjugate wave functions. The 
Hermitian conjugation of Hqq switches the left and right vectors. The ortho normality 
condition in the biorthogonal basis can be written as (^1^) = 8jk- 

The diagonalization of the effective SMEC Hamiltonian can be achieved by means of 
an orthogonal transformation. The resulting eigenvalues are real if all reaction channels 
are closed. Above the threshold, the transformation becomes non-unitary: 



and it yields complex eigenvalues. Physical resonances can be identified with narrow 
poles of the S"-matrix [121 E3 E2] , or using the Breit- Wigner approach, which leads to a 
fixed-point condition [20|[22|H3]. 

2.2. How to monitor effects of branch points and avoided level crossings in the complex 
energy plane ? 

The energy- dependent term (T2]) describing the coupling of closed quantum system 
eigenfunctions $j to the environment of scattering states does not act uniformly on 
all ISM states. A useful indicator of effects caused by an anti-Hermitian coupling term 
W I (E) on biorthogonal eigenfunctions {^k} of the effective Hamiltonian is the phase 
rigidity [H g5l 06] : 



given by the ratio of biorthogonal and Hermitian norms of an eigenf unction. The 
indicator (HJ) can also be written in a form j^6] : 



(3) 




(5) 
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which better illustrates its physical meaning. In this expression, the angle 6j arises from 
a transformation of so that Re'Q/j and Im^/j are orthogonal and it characterizes the 
degree to which the eigenfunction tyj is complex. 

Phase rigidity varies between 1 and 0. It equals 1 for bound state eigenfunctions. 
For unbound states, the condition rj=l means that the non-Hermitian coupling W J {E) 
exerts a negligible effect on the structure of an eigenfunction, i.e., its biorthogonal and 
Hermitian norms are identical. The energy variation of Tj is an internal property of a 
considered open quantum system and measures a mutual influence of the neighboring 
resonances. 
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Figure 2. The energy dependence of the phase rigidity for the wave functions of 
the Of and Of SMEC states of 16 Ne. The dashed curve shows the difference between 
energies of these states in the complex energy plane. 



Figure |2] shows the phase rigidity as a function of the scattering energy E for 
the two lowest-energy + levels of the effective Hamiltonian of 16 Ne. These SMEC 
calculations are performed in a pi/2, d 5 / 2 , S1/2 model space. For Hqq we use the ZBM 
Hamiltonian [17]. The residual coupling Hqp between the internal space and the 
surrounding continuum is generated by the contact force: Hq P = H P q = V 5(ri — r 2 ) 
with Vq = —1100 MeV fm 3 . The + ISM eigenstates are coupled through the common 
physical one-proton emission channels [ 15 F(I n ) <g) 1 H(lj)]^ l with I n = l/2 + ,5/2 + , and 
l/2~ which have thresholds at E — (elastic channel), 0.67 MeV, and 2.26 MeV, 
respectively. (Each physical channel, specified by a states I n and J n in daughter 
A ~ 1 X and parent A Y nuclei, respectively, may correspond to a number of different 
mathematical channels characterized by different (Zj)-values of a nucleon which are 
allowed by the coupling of I n and J 71 ".) These are all possible one-proton emission 
channels in 16 Ne, described in ZBM space. 

One can notice two local variations of the phase rigidity for Of and 0^ 
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eigenf unctions. The first one appears close to the elastic channel threshold whereas 
the second one, at higher energies, indicates the avoided level crossing between 0+ 
and O2" eigenvalues in the complex energy plane. Avoided level crossings are traces 
of exceptional points |HJ 09j |50] in the parameter space of the effective Hamiltonian 
[5T| 152"] . The dashed curve in Fig. [2] shows a distance in the complex energy plane 
between the two + eigenvalues. The minimum in Ae = \e\ — £2\{E) coincides with the 
minimum in r(E) at a position of the avoided crossing. However, neither maximum nor 
minimum of Ae coincides with the minimum of the phase rigidity around the elastic 
reaction threshold. 

A striking example of the configuration mixing caused by the coupling term W(E) 
is shown in Fig. |3j The phase rigidity for 2f (i = 1, • ■ ■ ,4) states in 16 Ne is plotted 
here for Vq = —1357 MeV fm 3 . The 2 + states are interconnected via the coupling to 
the inelastic channel: [ 15 F(5/2 7r ) (g) 1 H(/j)]^. One can clearly see that at this value of 
the continuum coupling strength there appears the degeneracy (exceptional point) of 
2g" and 2\ levels at which the phase rigidity drops abruptly to 0. The remaining 2 + 
eigenfunctions are mere spectators of this violent wave function rearrangement. At low 




1 — 1 — ■ — ■ — ■ — 1 — ■ — ■ — ■ — ■ — 1 — ■ — ■ — ■ — 1 — ' 

5 10 15 



E [MeV] 

Figure 3. The energy dependence of the phase rigidity for lowest 2+ (i = 1, • • • ,4) 
eigenfunctions of the SMEC Hamiltonian in 16 Ne. The 2^ and 2\ eigenvalues form an 
exceptional point at E ~ 4.8 MeV. 

energies, close to the branch point at the channel threshold (E=0.67 MeV), one can see 
a strong mixing of 2f and 2% eigenfunctions. 

2. 3. Near-threshold instability of the wave function 

Both the branch points associated with reaction thresholds and the avoided level 
crossings are essential elements of the configuration mixing in open quantum systems. 
Abrupt variations of the phase rigidity as a function of the energy for certain eigenstates 
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of the effective Hamiltonian are indicative of wave function instability. Another indicator 
of such instabilities could be the continuum-coupling energy correction, E corr (E) = 
($i\W(E)\$i), to the ISM eigenvalue ($i\HQQ\$i). FigureH (left) shows the continuum- 
coupling correction for the ground states of oxygen isotopes calculated for a single 
channel corresponding to a neutron coupled to the ground state in the daughter nucleus 
A ^0. 
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Figure 4. The continuum-coupling energy correction E co „(E). The value E=0 
corresponds to the first (one- neutron emission) threshold. For Hqq, we take USD 
Hamiltonian for the sd shell [S3] and KB' interaction for the pf shell [53]. The cross- 
shell interaction is the G- matrix |55j . The coupling to continuum is given by the 
Wigncr-Bartlett interaction with the strength Vq = —414 MeV fm 3 . Left: ground state 
configurations of 23_29 0. Only a contribution from a single channel corresponding to a 
neutron coupled to the ground state in the daughter ^ _1 nucleus is included. Right: 
the ground state configuration of 27 0. Contributions from couplings to various positive 
parity states in 26 are indicated. 



One can see that the continuum-coupling energy correction does depend on wave 
functions involved in the coupling matrix element, i.e., the wave function in the 
parent nucleus [N, Z] and the channel wave function [ z X A r_ 1 (/ ?r ) ® 1 \i{l . The 
correction E co „ varies strongly in a broad interval of energies around the neutron 
emission threshold: from ±4 MeV in 28 to ±2 MeV in 26 and 24 0. As a result 
of the centrifugal barrier, a maximum of |-E C orr| is shifted above the threshold if the 
angular momentum I of a neutron is different from zero. For 1=0, the first derivative 
of \E corr (E)\ is discontinuous at the neutron emission threshold. A similar shift can 
be caused by a Coulomb barrier if the continuum coupling involves charged particle 
emission channels. Another interesting observation is the odd-even staggering of E corr 
[40J: a blocking of the virtual scattering to the particle continuum by an odd nucleon 
diminishes the continuum correction to the binding energy in odd-iV nuclei. 

Figure H] (right) shows individual contributions to the continuum-coupling energy 
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correction in 27 O coming from couplings to various reaction channels [gOis^ 71 ") <8> 
1 n(/j)]^T. The largest contribution to E COTT comes from the coupling to the lowest 
physical inelastic channel: [gOis^^) ® 1 n(/j)] 3/ ' 2+ . The same wave appears in the elastic 
channel [ 8 0i 9 (3/2^) <g> 1 n(lj)} 0+ contribution to E corr in the neighboring nucleus 28 (see 
Fig. Hon the left). 

The angular momentum dependence and charge dependence of the continuum- 
coupling correction around the particle-emission threshold can be explained in the same 
way as a characteristic threshold behavior of scattering and reaction cross sections [57] 
(Wigner cusp). Due to the unitarity of the S- matrix and the resulting flux conservation, 
Wigner's threshold effect has a broad impact on various channel wave functions due to 
channel coupling [581 EH1 EH EH [62]. Hence, the branch point singularity at the opening 
of a given reaction channel can induce non-local correlations between wave functions in 
distant channels. This effect has been found experimentally, for example, in the coupling 
of the analogous channels in (d,p) and (d,n) reactions [63]. These channel-channel 
correlations are yet another manifestation of the eigenfunction correlations induced by 
the presence of branch points or exceptional points. 

Recently, it has been found that a threshold law, similar to the Wigner's law for 
scattering and reaction cross-sections, holds for one-nucleon overlap amplitudes [64], 
which probe shell occupancies; hence, they characterize many-body correlations. An 
immediate consequence of this finding is that the influence of the branch point on the 
configuration mixing in weakly bound mirror states will be different, even for identical 
separation energies in mirror states [65] . 

The size of the energy correction is a measure of the mixing of unperturbed ISM 
wave functions due to the coupling to the reaction channel(s), either open or closed. For 
ISM wave functions below all reaction thresholds, only the Hermitian part of the coupling 
term W(E) acts. For unbound ISM wave functions, the competition of Hermitian W R 
and anti-Hermitian W 1 parts is an essential ingredient of the formation mechanism of 
exceptional points in the spectra of resonances. As shown above, | -Ecorr I is both state- 
and system-dependent, i.e., the continuum coupling does not induce a global instability 
of all ISM eigenstates close to the particle emission threshold. On the contrary, for those 
ISM wave functions which are susceptible to the continuum coupling, the near-threshold 
instability is seen predominantly in a single pair of unperturbed ISM states. Other ISM 
states with the same quantum numbers remain spectators even though all of them are 
interconnected by the coupling to the same reaction channel(s). Both avoided level 
crossings (or exceptional points) and branch points associated with reaction thresholds 
induce the configuration mixing of continuum states. For weakly bound states, only 
branch point (s) can contribute to the configuration mixing. 

In general, the mixing mechanism associated with the coupling of ISM eigenstates 
to the reaction channel (s) enhances the similarity of a bound, near-threshold ISM state 
^(J 71 ") with the decay channel [ty^ £g> (r z lj)] J7r , i.e., increases the corresponding 
spectroscopic factor (\I/(J 7r )||a^||\I/ ( ^(P r )) 2 , where iff^ is the wave function of the 
daughter system. (The tilde symbol above bra vector signifies that the complex 
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conjugation in the dual space affects the angular part and leaves the radial part 
unchanged.) This 'alignment' of a near-threshold eigenstate of the effective Hamiltonian 
with the corresponding reaction channel has been demonstrated in the approximation of 
a one-nucleon continuum [56]. The phenomenon is based on three general principles: (i) 
the branch point nature of the reaction channel threshold, (ii) the unitarity of the 
S-matrix, and (iii) the 2x2 matrix structure in the eigenfunction mixing. These 
principles neither invoke any particular interaction in the internal space nor any special 
structure of eigenfunctions involved in the mutual coupling. Hence, the alignment of 
near-threshold eigenfunctions of the effective Hamiltonian with the reaction channel 
should be present in any open quantum many-body system which has complex reaction 
channels. In particular, the entanglement of ISM wave functions due to their coupling 
to a common particle emission channel should be responsible for the omnipresence 
of clustering phenomena in states close to their cluster decay thresholds. Indeed, as 
illustrated in Fig. [5j a-cluster states have been found among the ISM wave functions 
close to the a-particle emission thresholds in, for example, 12 C, whereas correlated two- 
neutron pair has been seen in the halo ground state of n Li close to the two-neutron 
emission threshold. 

a cluster states halo states 
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Figure 5. Cluster states close to decay thresholds. Left: Excited + state in 12 C 
(Hoyle state) just above the 8 Be+a threshold. Right: Subthreshold halo ground state 
of n Li. 



The threshold features are generic: they are neither accidental nor arise from 
particular properties of the nuclear Hamiltonian. This universality forms a foundation 
of Ikeda systematics jHS] which can be extended naturally to halo phenomena, see Fig. 
In this formulation, Borromean systems correspond to an example of clustering where 
the cluster itself is not a self-bound system. 

We are still far away from a quantitative understanding of the organization of 
ISM eigenfunctions near the particle emission threshold. The structure of the ground 
state wave function of 11 Li is a particularly interesting case-study due to an intricate 
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sequence of low-lying two- and one-neutron emission thresholds which are separated 
by less than 30 keV (!). This approximate degeneracy of two different emission 
thresholds may have important consequences on the structure of halo two-neutron pair. 
The relative weights of configurations, such as [ 9 Li(3/2~)] ® [n 2 (0 + )], corresponding 
to the two-neutron continuum of 9 Li, and {[ 9 Li(l/2~)] ® [n(l/2~)]} 1+ <g> [n(l/2~)], 
{[ 9 Li(l/2-)] ® K1/2+)]} 1 - ® [n(l/2+)], or {[ 9 Li(3/2~)] ® [n(l/2-)]} 2+ ® [n(l/2-)], 
associated with an unbound 10 Li, reflects the influence of these two almost degenerate 
branch points on the structure of the ground state wave function of 11 Li. 

The exceptional points, whose dramatic effects on the phase rigidity can be seen in 
Fig. [3l are also generic in Hamiltonian systems. They arise as a result of an interplay 
between Hermitian and non-Hermitian parts of the continuum coupling term W(E) for 
energies above the elastic channel threshold. The number of exceptional points depends 
only on some basic features of the system, such as its dimensionality and quantum 
integrability [67]. On the contrary, the precise position of any exceptional point in the 
spectrum of eigenvalues of the energy- dependent effective Hamiltonian depends on the 
choice of the ISM Hamiltonian Hqq, the strength of the coupling term W(E), and 
the energy E of the system. Hence, its finding provides the test of both the effective 
interaction and the configuration mixing in any OQS. 

It has been found recently that low-energy exceptional points appear for realistic 
values of coupling to the continuum and hence could be accessible experimentally 
[52] . At low excitation energies, they could be seen, for example, as individual 
peaks associated with a jump by 2tt of the elastic scattering phase shift. Also in 
the neighborhood of an exceptional point for avoided crossing of resonances, this 
characteristic imprint of the scattering phase shift remains, which gives a real chance 
that traces of the exceptional could actually be searched for in reaction studies. 

Complex and biorthogonal eigenstates of T-Lqq(E) provide a natural basis in which 
one can express the resonant part of any scattering wave function: 

t c s = ci + E«A- . (6) 

j 

where aj = ^4^^j^- and ^ = (1 + H P q)^> j. Dominant contributions to T% 
inside an interaction region are given by eigenfunctions tyj of Hqq, i.e., 

T C E ~E«^' • (7) 

j 

Both branch points (reaction thresholds) and double poles of the S matrix (exceptional 
points) lead to a non-separable entanglement of two eigenstates of the effective 
Hamiltonian which manifests itself in a singular behavior of matrix elements of any 
operator which does not commute with %qq [52], ETJ [68] • These unusual features are yet 
another facet of the profound change of the nature of ISM states close to the reaction 
thresholds and avoided resonance crossings which are responsible for clustering and 
strong mixing in OQSs. 
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2.4- Future developments: happy marriage of nuclear structure and nuclear reactions 

To understand the formation of a network of correlated states, e.g., the chain of isotopes 
from a well-bound 9 Li to an unbound 12 Li, a structure of the Hoyle resonance, and the 
radiative capture reactions, such 12 C(a,7) 16 0, a theoretical framework is required that 
would provide a unified approach to nuclear structure and reactions. Such a theoretical 
framework could be developed based on the GSM [231 EH ESI EE]- This model, formulated 
in the Rigged Hilbert Space j28l EH] and using a complete Berggren ensemble [27], is 
a generalization of the nuclear ISM for a description of bound states, resonances, and 
many-body scattering continuum. In a nuclear structure application, solutions of the 
GSM can be found by diagonalizing a complex-symmetric Hamiltonian matrix. The 
'dimensional catastrophe' in GSM when increasing the number of active particles is 
much more serious than in the standard ISM because each single-particle continuum 
state in the Berggren ensemble becomes a new shell in the many-body GSM formulation. 
This acute problem has been alleviated by recent progress in the generalization of the 
Density Matrix Renormalization Group [EH [TO, E] method to non- Hermit ian, complex- 
symmetric matrix problems [721 E2] . 

Significant progress has also been made in applications of realistic interactions in 
GSM [71]. Finally, powerful techniques for a selection of physical resonances, based on 
the overlap method, have been developed [231 EE] • Still, much work is needed to develop 
appropriate effective interactions which would allow a systematic study of the structure 
of nuclear states in long isotopic or isotonic chains: from the valley of stability towards 
particle drip lines. 

So far, most applications of the GSM have addressed nuclear structure phenomena. 
Further progress can only be achieved if the method could be fully extended to reaction 
problems. This can be achieved by the coupled-channel formulation of the scattering 
problem: 

/ *-XE< c ; r l*»> l(*n\H - E8 nm \y m )] (^rn\c';r'}Q c/ (r') = . (8) 

q> m,n 

A similar approach, based on the Resonating Group Method [7_5l [7_6j [77], has been 
developed in the no-core shell- model approach to nuclear scattering problems [7_8j [79] . 
In the above equation, H nm = (^> n \H\^> m ) is the GSM Hamiltonian matrix and 
O c - n = (c;r\^f n ) are overlaps of the reaction channel and the GSM many-body state. 
Formally, Eq. (jSJ) takes the same form independently of the choice of reaction channels; 
the whole complication is hidden in the overlaps C]n . Their calculation is a formidable 
task for the GSM approach. It requires a complicated change of the coordinate system 
which makes the calculation difficult for heavier projectiles, even for reactions with two- 
body asymptotic conditions. Any future progress in the application of GSM to nuclear 
reactions is ultimately related to the progress in the development of new algorithms for 
fast calculation of overlaps. 
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